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We in the paper study the metric perturbations generated in a bouncing universe driven by the 
Quintom matter. Firstly, we review the background evolution of Quintom Bounce and the power 
spectrum of scalar perturbations. Secondly, we study the non-Gaussianity of curvature perturbations 
and then calculate the tensor perturbations of the model. 
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I. INTRODUCTION 

Inflation was invented to resolve problems existing in 
hot Big Bang cosmology, such as flatness, horizon, pri- 
mordial monopole problem[l[ (for some early attempts 
see Refs. Q). After more than twenty years' devel- 
opment, one has obtained a deep sight at this theory. 
However, it is still puzzled by an initial singularity which 
exists in usual inflationary models Q. 

It has been suggested that, bouncing cosmology, which 
requires our universe initially experience a contracting 
phase before the hot Big Bang expansion, could provide 
a solution to the problem of the initial singularity. For 
years, models of a bouncing universe have received a lot 
of attentions and there have been a number of works on 
constructing this scenario. For example, there are mod- 
els with singular bounce such as the Pre-Big-Bang Q 
and cyclic/Ekpyrotic p|; also in Refs. [1, 0> @] some 
non-singular bounce models were considered where the 
gravitation action was modified by higher order correc- 
tions; and, see Ref. for a recent review on various 
models of bouncing cosmology. 

Recently, a new bounce model has been proposed [l(| , 
dubbed Quintom Bounce. In this model, a bouncing uni- 
verse was obtained within the standard 4-dimensional 
Friedmann-Robertson- Walker (FRW) framework by 
making use of Quintom matter [111]. The key feature of 
this model is that, the null energy condition (NEC) has 
to be violated for a short while around the bounce point 
and after that the equation-of-state (EoS) of our model 
w is able to transit from below —1 to above —1 which 
makes the universe enter into normal expanding history. 

The model of Quintom Bounce has presented a very 
interesting picture of the early universe. Firstly, the 
investigation of the cosmological perturbations in this 
model [l2] has shown a combination of some aspects 
found in some recently studied non-singular bounce mod- 
els II 3. 1 141 and some others in singular bounce models 
(HI. Il6l TlTl [l8j . Secondly, in a recent study on Quin- 
tom Bounce, the authors of Ref. [li| have found that a 
concrete model of Quintom Bounce are able to provide a 
scale-invariant spectrum in ultra-violet regime, and also 
give rise to an oscillation signature which could be veri- 
fied in the forthcoming astronomical observations. 

One important lesson of studying the primordial curva- 
ture perturbations is to investigate its bispectrum which 



describes the non-Gaussianity of the power spectrum [20j]. 
There is a number of information coded in this bispec- 
trum, such as magnitude, shape, sign, and even running. 
To probe these signatures, we expect to distinguish var- 
ious models of the very early universe. For example, in 
the usual slow-roll inflationary model, it is pointed out 
that non-Gaussianity is negligible due to the suppression 
of slow-roll parameters 12 l l . 22 )1 ; however, in the models 
of Ekyrotic/cyclic [23], [2J, [25| or island universe [2rjj]. a 
large non-Gaussianity is predicted. Observationally, cur- 
rent cosmological data(2?], [28| is consistent with Gaus- 
sian distribution, however, the forthcoming observations 
will be sensitive to the non-Gaussianity with much higher 
precision, such as Planck satellite [29(. Given the consid- 
erations above, we in this paper calculate the bispectrum 
of a Quintom Bounce model. Our results show that non- 
Gaussianity in this model is still suppressed by slow-roll 
parameters, but there is an interesting oscillation signal 
on the non-linear parameter and its maximal value is 
mildly bigger than that in the usual scenario of slow-roll 
inflationary models. 

Another important clue to discover the information of 
the very early universe is to study the relic gravitational 
wave background (GWB) formed by primordial tensor 
fluctuations. There have been a number of detectors op- 
erating for the signals of primordial GWB, e.g. Planck 
[U, Big Bang Observer (BBO) LIGO 0|. More- 
over, the indirect detection of GWB attracts a lot of in- 
terests of the next gen eration of CMB observations, see 
related analyses [32ll33l.r34l.r35l|. The basic mechanism for 
the generation of primordial GWB in cosmology has been 
discussed in Refs. [36], [37j]- Usually, inflation theory pre- 
dicts that the power spectrum of primordial tensor fluc- 
tuations is scale-invariant and its value is roughly equal 
to scalar spectrum times a slow-roll parameter defined as 
e = — H/H 2 ]3§i. |39|. However, this so-called consistency 
relation is not neces sary to be valid in a bounce scenario. 
For example, Ref. [4f| has investigated the primordial 
gravitational waves in singular bounce models and pre- 
dicted an undetectable GWB. We in this paper study 
the relic gravitational waves in Quintom Bounce with 
Coleman- Weinberg potential. Interestingly, we find that 
the power spectrum of primordial gravitational waves is 
scale-invariant both in ultra-violet and infrared regimes 
but strongly oscillate in the middle band which is related 
to the bounce, and also find a sunken signature on the 
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energy spectrum. These results on one side support the 
conclusions of Ref. in which the scalar perturbations 
were discussed, and on the other side provide another 
approaches to searching the signals of a bounce. 

The outline of this paper is as follows. In Section 
II, we briefly review the basic scenario of a Quintom 
Bounce model with a the Coleman- Weinberg potential, 
and the curvature perturbation. In Section III, we make 
the calculation of non-Gaussianity in the model of Quin- 
tom Bounce, and find that the non-linear parameter /jyx 
is oscillating around a small value which coincide with 
that in usual inflation model. In Section IV, we inves- 
tigate the tensor part of gravitational perturbations in 
this model, and then give the power spectrum and spec- 
tral index correspondingly. Considering the influences of 
possible factors during the late time evolution, we discuss 
the transfer function for the gravitational waves and fi- 
nally obtain the energy spectrum of GWB we may probe 
in future observations. Section V contains discussion and 
conclusions. 



II. A REVIEW: BACKGROUND AND LINEAR 
PERTURBATIONS 

As proposed in Ref. [1(3] , a model of Quintom Bounce 
can be described in terms of scalar fields which minimally 
couple to the four dimensional Einstein's gravity. The 
simplest model is realized by two scalar fields, with its 
Lagrangian given by 

C = -\d„ <pd»<p + \d^d^ - V{<f>, i>) , (1) 

in a spatially flat Friedmann-Robertson- Walker (FRW) 
universe where the signature of the metric is taken to be 
(— , +, +, +). Here the essential component is the scalar 
field ip. It plays a crucial role in giving a bouncing so- 
lution smoothly. Without it, the model in (1) will be 
similar to the traditional inflation with a single scalar 
field, which as we know suffers from the problem of the 
initial singularity. 

In collaboration with Qiu, Xia and Li[l^|, we take the 
potential to be only the function of the field <f>, and specif- 
ically of Coleman- Weinberg form [41] : 

which takes its maximum value Xv /16 at <p = and van- 
ishes at the minima when <p — ±v. Therefore, the scalar 
field ip merely affects the evolution around the bounce 
but decays out quickly when away from the bouncing 
point. 

This model performs an interesting evolution for the 
background universe [19]. Initially, (p stays at its left vac- 
uum — v and the kinetic term of ip is sufficiently small. 
This looks very natural and there is no argument for any 
fields being outside the Planck scale, so it does not suffer 



the initial condition problem which appears in the usual 
inflationary model [42j,|43( . At the beginning of the evolu- 
tion, the field <p oscillates around the vacuum point — v, 
so the EoS of the universe oscillates about w = and 
averagely the state looks similar to a matter dominated 
one. Since the universe is contracting, the amplitude of 
the canonical field oscillation becomes larger and larger, 
in the meanwhile the contribution of ip also becomes im- 
portant. When the field reaches the plateau, the energy 
density of (p would be cancelled by that of ip and so the 
bounce happens. After the bounce, as the field (p moves 
forward slowly along the plateau, the universe enters into 
a slow-roll phase and the EoS of the universe is approxi- 
mately — 1, very much alike inflation. The process to link 
the contraction and expansion is a smooth bounce, dur- 
ing which the evolution of the hubble parameter can be 
treated as a linear function of the cosmic time t approx- 
imately. Finally, when the field <p "drops" into the right 
vacuum +v, it will oscillate again and reheat the uni- 
verse. To present the above analysis clearly, we give the 
numerical calculation of background evolution in Fig. [TJ 
From this figure, we can see that, although the potential 
is symmetric w.r.t. the field <p, the background evolution 
is asymmetric w.r.t. the cosmic time. 
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FIG. 1: A plot of the evolution of the EoS w, the hubble 
parameter H and the canonical field 4> in the model ([TJ where 
a Coleman- Weinberg potential is considered. In the numerical 
calculation we choose the parameters A = 8.0 x 10~ 14 , v = 
0.82Af p ;, and the initial conditions as: <p = — 0.82M P ;, <p = 
3.0 x 10- 10 M^, ip = -0.72M pi , ip = 5.0 x 10" 13 M^ where 
M pl = 1/VG. 

Note that, the asymmetry of the evolution is a sig- 
nificant character since we are able to obtain a scale- 
invariant spectrum in virtue of this asymmetry. The rea- 
son is as follows. Since the evolution of the universe is 
asymmetric w.r.t. the bounce point, it is possible for the 
primordial fluctuations keeping inside the horizon when 



3 



the universe is contracting. And if so, while these per- 
turbations are able to escape to the super-hubble region 
during the inflationary period, the initial Bunch-Davies 
condition for them can be basically saved and transferred 
through the bounce. Therefore, this scenario provides 
a possible approach to obtaining a scale-invariant spec- 
trum. However, as will be shown in the following, we 
find that the sub-hubble perturbations still deviate a lit- 
tle from the pure incoming plane wave on the matching 
surface between the bouncing phase and the expanding 
phase. This deviation would bring some wiggles on the 
corresponding scale of the primordial power spectrum, 
which has been shown in Ref. [HI can be detected by 
the future observations. 

Now we take a brief review of the linear curvature per- 
turbation (for details we refer to Ref. [19j |. and see Ref. 
(4^ | for a comprehensive review), then study its evolu- 
tion in the model of Quintom Bounce and obtain the 
expression for the primordial power spectrum. Under 
the longitudinal (conformal-Newtonian) gauge, the met- 
ric perturbations are given by 

ds 2 =a 2 (T)[~(l + 2$)dT 2 + (l-2^)dx t dx i ] , (3) 

where we introduce the comoving time r defined by 
dr = dt/a. We start with the equation of motion of 
the gravitational potential 

$" + 2(H - %-W + 2{H' ~ U^--)^ - V 2 $ 

= 8ttG(2H + ^-tWH , (4) 

9 

where TL = a' /a is the comoving hubble parameter and 
the prime denotes the derivative with respect to the co- 
moving time. This equation can be derived from the ba- 
sic perturbation equations directly (we refer the complete 
derivation to Ref. [Hj]). 

As we have pointed out in describing the background 
evolution, the energy density of the scalar ip is usually 
negligible away from the bouncing point and hence we 
have ip' ~ 0. Therefore it is a good approximation to 
neglect the right hand side of Eq. (01 when the universe 
is far away from the bounce. Moreover, we find that 
there is the relation 2H. + 4t- ~ in the bouncing phase 
and so the perturbation of ijj decouples from Eq. (@|, 
which has been checked in Ref. [l2j]. Thus, in the paper 
we will neglect the r.h.s. of Eq. (g]), and just focus on 
the adiabatic fluctuations which can be determined by a 
single physical scalar field degree of freedom $. 

If the evolution of the background is known, all other 
perturbation variables can be determined from <£>. A fre- 
quently used variable is the curvature perturbation in 
comoving coordinates C which is defined as, 
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7 ($' + ft$) , 



(5) 



and this variable can be calculated from $ and the back- 
ground parameters. For example, when the universe is in 



a nearly de-Sitter like expansion, there is a simple rela- 
tion £ ~ $/e with the slow-roll parameter e = —H/H 2 . 
In usual case this variable is well known to describe the 
adiabatic perturbations on large scales since it is a con- 
served quantity on super-hubble scales according to the 
equation 



9(1 + w)H 2 



(6) 



and its dynamics can be simply and conveniently de- 
scribed by the equation of motion for Mukhanov-Sasaki 
variable [45|]. However, this equation becomes ill-defined 
both at the bounce point with a vanishing hubble pa- 
rameter H = and the cosmological constant bound- 
ary w = —1. This point has been remarked in Ref. 
O EH Ell • Consequently, we investigate the evolution 
of the gravitational potential $ directly in deriving the 
curvature perturbation, and then moves to £ after the 
universe entering the expanding phase. 

As introduced in the above, the universe in this model 
experiences three phases, a contracting one, a bounce, 
and finally an inflationary one. We can resolve the per- 
turbation equation in each phase respectively. In the de- 
tailed calculation, we take the Bunch-Davies vacuum as 
the initial condition $k ~ ^|== \4>\e~ lkr , and then make 
those solutions smoothly pass through the linking point 
applying the matching conditions [47| (see also [IR l49l. [Hcl] 
for a recent study) . The detailed derivation has been pre- 
sented in Ref. [19j |. and here we give the dominant part 
of the final gravitational potential directly: 



AirG 



-ikr 



1 - 



m B -e lk,HB ^ 



2k 



sm[k/Ti.B^ 



(7) 



where 7is_ and TLb+ represent the comoving hubble pa- 
rameter at the beginning and the end of the bouncing 
phase respectively. By comparing the coefficient ([7]) and 
the initial form of obviously the dominant mode of the 
gravitational potential deviate from the Bunch-Davies 
initial condition when the inflationary phase takes place. 
Moreover, we have the approximate relation Q ~ $/e. 
So we eventually have the expression of the primordial 
power spectrum for the curvature perturbation 1 
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(8) 



From this result, one obviously find that the first term 
provides a nearly scale-invariant spectrum which is con- 
sistent with current cosmological observations. However, 



1 This result seems similar to an inflation theory with its initial 
condition as a- vacuum |5 ll. [53| . However, these two are different 
in the detailed characteristics of the perturbations and the basic 
mechanisms for the generation of those perturbations. 
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the second term shows that there is a wiggle on the spec- 
trum, due to the modified initial condition by the bounce 
relative to the standard inflation. Apparently, this os- 
cillation term could affect the CMB temperature power 
spectrum and LSS matter power spectrum [l9j|. 



III. BISPECTRUM AND THE NON-LINEAR 
PARAMETER f NL 

In this section we investigate the evolution of the non- 
linear part of scalar perturbations in our model. As we 
have discussed in the last section, the universe experi- 
ences a nearly exponential expansion after the bounce, 
and so we can deal with the curvature perturbation of 
Quintom Bounce as in inflation theory except for the 
initial condition being modified. Thus, the slow-roll pa- 
rameters can be defined in our model, with e = —H/H 2 

To expand the action to the third order of £ and drop 
the terms suppressed by slow- roll parameters e and 77, we 
have the final cubic action for the curvature perturbation 
during the inflationary phase plj. 



S, 



dtd 3 x 



ga 2 t 2 d- 2 t + 2/(C)f li 



where 



/(C) 



(9) 



(10) 



Note that the last term ^ 1 1 in the third-order action 



and g = 4a 3 e 2 H. 
Note that the 

can be absorbed by field redefinitions of (. It is taken as 

C = C + /(0- (11 



This field redefinitions do not affect any other of the 
0(C 3 ) terms in the third-order action, since the term / 
are quadratic in Moreover, the last term in the third- 
order action can be cancelled by one extra quadratic part 
of this field redefinitions which is proportional to the first- 
order equations of motion ^ 1 1 exactly. 

After deriving out the third-order action, now we are 
able to get the three point function for our model. It 



can be computed using the path integral formalism in 
the interaction picture as follows 

<C(t,ki)C(t ) k 2 )C(t,k 3 ) >= 

ij^ dt< [C(t,ki)C(«,k 2 )C(i ) k 3 ),i 3 («)] >, (12) 

where U denotes an initial time for the modes deep in- 
side the horizon, and L3 is the third order perturbative 
lagrangian. 

Recall in Eq. ([5]) we have obtained the linear curvature 
perturbation which deviates from that in the standard 
inflation theory by multiplying an extra factor defined as 
follows, 



1 3Hb - ■ 2k 
C (k) = 1 — — sin ; 
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(13) 



Considering this factor in the calculations, the leading 
order contributions to the three point correlator in our 
model are listed in the following, 

(a): contribution from ( 2 d~ 2 (. To make the expression 
simplified, we define — §G 2 | which is just the curva- 
ture spectrum in inflation models and K = ^ . ki . So we 
have the three point correlator from £ 2 d~ 2 (: 



(2 7 r)^)(^k l )(P c ) 2 (n^) x -^£^ 2 • (14) 
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(b): contribution from the field redefinition. 
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i C{ki) 8 ^i^j C(ki) 



(15) 



Moreover, since non-Gaussianity measures the devia- 
tion of CMB power spectrum from the Gaussian distribu- 
tion, we can define a non- linear parameter f^L as follows, 



c = c 9 + ^mc;- <c;>) ■ 



(16) 



So we eventually have /jvl to characterize the size of 
non-Gaussianity, 



/jVL — -TT 
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(17) 



There are two limiting cases of non-Gaussianity, which ^ » ij). To the case of equilateral form, /jvx is given 
are of particular interests for observations. These are 
equilateral form (k\ ~ ki ~ fca) and local form (fci ~ 
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IV. GRAVITATIONAL WAVE BACKGROUND 
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C(k) 



(18) 



where we have taken the limit k — k\ — &2 = k%. For 
the local form, which corresponds to that k% mode exits 
horizon much earlier than the other two, we have the 
non-linear parameter 



rlocal 
J NL 



2 + C(fc) e-r? 



(19) 



where the limit is taken as k = k\ = &2 3* One may 
notice that, the above results reduce to the single scalar 
slow- roll inflationary model when C(k) — > 1. However, 
since the factor C(k) has modified the initial condition 
of curvature perturbation when the universe enters infla- 
tionary stage, it can bring an oscillation signature on the 
size of non-Gaussianity as well. In order to compare our 
result with the non-Gaussianity predicted by usual infla- 
tion model, we plot Jnl of equilateral and local forms in 

Fig.m 




FIG. 2: The plots of the non-linear parameters /ml with 
different forms predicted by our model are shown by the ma- 
genta solid lines. To make the comparison, we also give Jnl 
predicted by a single scalar inflation model in blue dotted 
lines. 

Our results show that the non-Gaussianity in the Quin- 
tom Bounce model is still suppressed by the slow-roll pa- 
rameters. However, there is an oscillation signature on 
Jnl and the maximal value of /jvl is bigger than that 
in single scalar slow-roll inflationary models. The reason 
for this effect is that the dominant modes of the curva- 
ture perturbations have deviated from the Bunch-Davies 
form when they pass through the bounce and enter the 
inflationary stage. This is similar to the cases with non- 
Gaussianity generated from a modified initial condition, 
for example see Refs. [H, IBH [55lj . 



Now we turn to consider the evolution of gravitational 
wave background from the tensor part of the primordial 
metric perturbations. In order to standardize the deriva- 
tion, we use the same convention as in Ref. [5y|. To 
begin with, we give the metric containing the tensor per- 
turbations in the flat FRW background as follows, 

ds 2 = a(r) 2 [-dr 2 + {6ij + h ij )dx i dx j , (20) 

where the Latin indexes represent spatial coordinates. 
Here the tensor perturbation hij satisfies the following 
constraints: 



(21) 



Due to these constraints, we only have two degrees of 
freedom in hij which correspond to two polarizations of 
gravitational waves. 

By adding the anisotropic part of the stress tensor <7jj , 
we have the equation of motion for tensor perturbations, 



2^ 
a J 



V 2 hij = 16nGa 2 <Jij . 



(22) 



The Fourier transformations of the tensor perturbations 
and anisotropic stress tensor are give by, 

hij (r, x) = VT^G / - —j- Hij (r, k)e ikx , (23) 



(2tt)4 
d 3 k 



/d A k 
-£- T E ii (r,k)e ikx 
(27T) 2 



(24) 



Note that, what we are usually interested in are the 
distribution of the spectra of gravitational waves and the 
corresponding spectral index. Based on the above for- 
malism, the tensor power spectrum can be written as, 



P T {k,T) 



d(0\h%\0) 
d In fe 



327rG 7 



W l^(r,k)f , (25) 
and the definition of tensor spectral index n-r is given by 

nT = ^k ■ (26) 

The GWB we observed today is characterized by the en- 
ergy spectrum, 

1 d(0\ PGW (r)\0) 



&Gw{k,T) 



Pc(r) 



din A: 



(27) 



where pgw( t ) is the energy density of gravitational 
waves, and the parameter p c (r) is the critical density of 
the universe. In respect that the GWB we observed has 
already reentered the horizon, the modes should oscillate 
in the form of a sinusoidal function. Consequently, we can 
make use of the Friedmann equation H 2 (r) = ^^Pc(t) 
and then deduce the relation between the power spec- 
trum and the energy spectrum as follows, 

^(fc,r)^l_^_p T(fc;r)j (28) 
which will be used in the following calculations. 
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A. Tensor perturbations 

Now we follow one Fourier mode of the tensor pertur- 
bations, labelled by its comoving wave number k, and 
find that there are two paths which are different in the 
times of crossing the hubble radius. The evolution of 
tensor perturbations is sketched in Fig. [31 Initially all 
the perturbations stay inside the hubble radius in the 
far past. Since the hubble radius shrinks in the con- 
tracting phase, those modes with small comoving wave 
number exit the hubble radius while the large k scales 
still keep inside. When the bounce takes place, all the 
perturbations will keep inside the hubble radius because 
at that moment the hubble radius diverges. After that 
the bounce is followed by a slow-roll expanding phase, so 
these Fourier modes will escape out if the efolds for the 
post-bounce inflationary period is large enough. After 
that, these modes will reenter the hubble radius at late 
times after the slow-roll expanding phase has finished. 




Large k Small k 



FIG. 3: A sketch plot of the evolution of tensor perturbations 
with different comoving wave number k in our model. 

Therefore, we can classify the tensor perturbations 
with different comoving wave numbers k to two cate- 
gories, as the two lines sketched in Fig. [3] The blue line 
denotes a mode with a scale k which is large enough to 
keep it inside the hubble radius through the contracting 
phase and the bounce, and then escape outside during 
the post-bounce inflationary period; the green line con- 
sists of a mode with small k so that it exits the hubble 
radius in the contracting phase, then is pushed inside 
during the bounce and soon is pulled outside before the 
slow-roll expanding phase happens. 



Due to the symmetry of , we can express the two 
polarizations as one function v = ~J^{Hn + iH^-i). Ne- 
glecting the anisotropic stress tensor in the very early 
universe, we obtain the equation of motion for v: 

a" 

v" + k 2 v v = . (29) 

a 

Following the background evolution of the universe, we 
obtain three solutions of gravitational waves similar to 
what we did with scalar perturbations. For the universe 
which is contracting with its EoS oscillating around w = 
0, we have 

v = (T-f B .)i(AlH^[k(T-r B .)} 

+^4 2) [fe(r-f B _)]| , (30) 

where tb- = tb- + 2/Hb-- Here and are the 
p-th Hankel function of the first kind and second kind 
respectively. Besides, the parameters A\ and can 
be determined by the initial condition for gravitational 
waves, which is usually taken as Bunch-Davies vacuum 
v ~ e~ ikT /V2k. So we have A\ — and B% = -y^F/2. 
Therefore, the asymptotic forms of the solution to the 
tensor perturbation in the contracting phase is 

!3 
— iJ£j=(j- — tb-) 2 , outside horizon; 
(31) 
_i e -*k(T— tb-) i ns ide horizon. 

When the universe undergoes the bouncing phase, we 
have the approximate relation that — ~ ^aa 2 B = |. To 
solve Eq. (j2"9"j) . we have 

v(k, t) = 

( Cl cos[Z(r - t b )] + Dl sin[Z(r - t b % k 2 > | ; 

{ (32) 
{ Cle l ^- T B>+Z>jT e -J(T-T b) 7 k 2 < i j 

where we define I 2 = \k 2 — ||. Since the hubble parame- 
ter approaches zero when the universe is bouncing from 
a contraction to an expanding phase, all the modes of 
the perturbations would return to the sub-hubble region. 
However, from the above solution we interestingly find 
that, k 2 h (~ k 2 /a 2 B ) and H(~ a) are comparable. 

After the bounce, the slow-roll expanding phase takes 
place which drives the universe to inflate like a de-Sitter 
spacetime. In this case, the solution to the gravitational 
waves is given by 

V = (T-TB+)* >< \ElH^[k{T-T B+ )] 

+F£HW[k(T-T B+ )]\ , (33) 
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where v = ~ 
form, 



Y^— — | • This solution has an asymptotic 



k-i{T-f B +)- 1 {E t 



Ft), 



(34) 



after the modes exit the horizon. 

Having obtained the solutions of the tensor perturba- 
tions in different phases, now we need to match these so- 
lutions and determine the coefficients Cj, , E% and 
respectively. This procedure is much similar to the 
matching process of scalar perturbations as done in the 
previous section. For a non-singular bounce scenario such 
as the Quintom Bounce model, the continuity of back- 
ground evolution implies that both v and v' are able to 
pass through the bounce smoothly. So we match v and 
v' in (|3"Tj) and ([3"!?]) on the surface tb- , and those in (13^)) 
and (|33[) on the surface tb+- With these matching con- 
ditions, we can determine all the coefficients and finally 

However, as what we have analyzed at the beginning of 
this section, there are two paths for the tensor perturba- 
tions to evolve from a contracting phase to an expanding 
phase. So there are two possible results for E% — F^ . For 
the first case, the comoving wave number is large enough 
so that the tensor perturbations have never escape out- 
side the hubble radius, thus we have 



\El 



where a 



n a 1 
1 e 

1 4 
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2ik8rB 



)l , (35) 



sga-, and 5t b 



2p-, auu u, B — t b + - t~b— . For the second 
case where the modes of gravitational waves are in small 
k region, the expression E/£ — is given by 
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FIG. 4: The red solid curve represents the primordial power 
spectrum Pt and the spectral index nr of tensor perturba- 
tions in our model. The blue dash curve give the primordial 
power spectrum and the spectral index in a single scalar in- 
flation model. In the figure, we take the values of parameters 
the same as in Fig. Q] 



caused by the modified dispersion relation of the ten- 
sor perturbations when they pass through the bouncing 
phase 2 . However, when the comoving wave number k 
gets even smaller, the power spectrum is able to climb 
up and finally reaches a certain value with its spectral 
index returning to zero again. 



E 



F. 



y^ n B -{2l+HB-) - lSTB 

6 

8 ~~ 7^.B+) 



(36) 



Based on the above analysis, now we are able to derive 
the primordial power spectrum of gravitational waves. 
From the definition of Eq. (|25p . the primordial power 
spectrum is given by 



P T {k) 



UGH 2 



El 



(37) 



From Eqs. (f3"5| and (f3"6"|) , we can read that the spectrum 
are scale-invariant both at the large k and small k region, 
but oscillate when k is near to a critical value \pk- To 
illustrate the above analysis clearly, we do the numeri- 
cal calculation and plot the results of primordial tensor 
power spectrum and the corresponding spectral index in 

Fig. m 

One can see that in Fig. [4] when the value of k is large 
enough, the red solid lines converge at the blue dash lines 
with an oscillation. The amplitude of the oscillation gets 
the largest value when k approaches the neighborhood 
of the critical value y\, and soon drop down to a min- 
imal value when k gets smaller. This damping effect is 



B. Energy Spectrum of Today's GWB 

In the above section we discussed the behavior of ten- 
sor perturbations exhibited in primordial power spec- 
trum and the spectral index. However, we are more 
interested in how to recognize these perturbations in 
the GWB nowadays. Since the primordial gravitational 
waves are distributed in every frequency, once the effec- 
tive co-moving wave number is less than aH, the cor- 
responding mode of gravitational waves would escape 
the horizon and be frozen until it reenters the horizon. 
The relation between the time when tensor perturba- 
tions exit the horizon and the time when they return 
is a ou tH out = cii n Hi n . Therefore, we have the conclusion 
that, the earlier the perturbations escape the horizon, 
the later they re-enter it. Moreover, once the effective 
co-moving wave number is larger than aH , the perturba- 
tions begin to oscillate like the plane wave, as shown in 



2 A similar scenario of the primordial gravitational waves has been 
considered in Ref. |56H . where the authors have considered the 
damping effects from the spacetime noncommutativity. 
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Fig. El 

To relate the power spectrum observed today to the 
primordial one, one can define a transfer function T(fc, r), 
given by Refs. [H,[I3: 



T(k,r) 



where a 



0.80313 (1 + z{t) 



2tt VI 



Zk 



T 2 (a 



\)(-f a 
I a 



(38) 



i ; Q is determined by the EoS of the universe, 
z(t) is the redshift at the moment r and Zk is the redshift 
when the k mode of gravitational wave reenters the hori- 
zon. Here the factor 0.80313 comes from the damping ef- 
fect of freely streaming neutrinos [58|]. Moreover, the fac- 
tor ( ) 2 describes the redshift-suppressing effect on 
the primordial gravitational waves. The rest factor shows 
that, when the gravitational waves reenter the horizon, 
there is a "wall" lying on the horizon which affects the 
tensor spectrum. 

Considering today our universe is dominated by dark 
energy of which the EoS is w ~ — 1, we are able to obtain 
today's transfer function. Then we can get today's tensor 
power spectrum 



P T (k,T ) = 



0.80313 64GH? 

2(1 + Z k ) 2 7T 2 



El. 



(39) 



where Hi represents the hubble parameter in the infla- 
tionary stage. Eventually, the present energy spectrum 
of GWB is given by Qcwik, tq) — (aoH ) 2 ^rik, To). 

In Fig. [5] we plot the numerical results of the energy 
spectrum in our model. One can see that, when the fre- 
quency of GWB is large enough, the tensor energy spec- 
trum of our model would agree with the prediction of the 
single scalar inflationary model. However, when the fre- 
quency becomes smaller, the physics of a bounce begins 
to affect the behavior of the GWB. Moreover, there is an 
interesting sunken area in the middle band. This sunken 
signal is resulted from the modified dispersion relation 
of the tensor perturbations when they pass through the 
bouncing phase. 



V. DISCUSSION AND CONCLUSIONS 

Bouncing cosmology, due to the avoidance of the initial 
singularity, has attracted a lot of interests in the litera- 
ture [H M, EH, HI, [H, @ El (and see Ref. @ for a 
recent review). However, since it happens in extremely 
high energy regime, we hardly observe a bounce by ex- 
periments directly. So it is a debate whether a bounce 
has taken place or not. To find the evidences of a bounce, 
we need to know what can a bounce leave for observa- 
tions. This question is still discussed drastically in the 
literature, and one potential clue is to study the pri- 
mordial gravitational fluctuations. In the context of the 
Pre-Big-Bang scenario and in the cyclic/Ekpyrotic cos- 
mology, the primordial curvature perturbation strongly 
depends on the physics at the epoch of thermalization, 



«GW h 



10""' 
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10" ,s 
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10" 
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Energy spectrum of GWB - 




in Quintom Bounce 




Energy spectrum of GWB 1 




in inflation 







FIG. 5: The red solid curve represents the energy spectrum 
Qqw in our model. The blue dash curve gives the energy 
spectrum in a single scalar inflation model. In the figure, we 
take the values of parameters the same as in Fig. [T] 



and thus an uncertainty of a thermalized surface is in- 
volved [HI, [H, [13, d, HH ■ In the frame of loop quantum 
cosmology, it is argued that fluctuations before and af- 
ter the bounce are largely independent [671 ] (yet see Ref. 
[68| for some criticisms). We in this paper have stud- 
ied the perturbation theory of a Quintom Bounce model 
detailedly and show that there are some imprints of the 
bounce on CMB observations at large scales. In the main 
content we have analyzed both the linear and non-linear 
evolutions of scalar modes, and the tensor perturbations 
are also considered. 

The model we considered is constructed by a double 
field Quintom model with a Coleman- Weinberg poten- 
tial. We firstly have reviewed the background dynam- 
ics of this model, and obtained a scale-invariant scalar 
spectrum in virtue of an asymmetry of the background 
evolution around the bounce point. A similar but more 
phenomenological scenario has been studied in Refs. [69j] 
as a possible solution to the suppressed low multi-poles of 
the CMB anisotropies. Moreover, since the gravitational 
perturbations in sub-hubble region would change their 
propagations when pass through the surface between the 
contracting phase and the bounce, there would be an os- 
cillation signature generated both on the linear scalar 
modes and non- linear ones. We have also calculated 
the non-Gaussianity and shown that the maximal value 
of non-linear parameter /jvz, predicted by our model is 
mildly bigger than the usual one in single scalar slow- 
roll inflation, but the central value is still suppressed by 
the slow-roll parameters. So we expect that a large non- 
Gaussianity might be generated by other mechanisms [70l] 
in the frame of Quintom Bounce. 

We in the last part of this paper focus on the behavior 
of the gravitational waves in Quintom Bounce. Due to 
the effects of a bounce, the solution of the tensor per- 
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turbation is quite different from the usual one. In our 
analysis, we find that the physics of a bounce would af- 
fect the evolution of primordial tensor perturbations at 
large scales of the universe, which corresponds to the 
physics in very early time. The behavior of the energy 
spectrum of the GWB in our model is similar to that 
in the single scalar inflationary model in high-frequency 
regime. However, for low-frequency regime the differ- 
ence becomes larger. Moreover, there is a sunken area in 
the middle band which links high-frequency regime and 
low-frequency regime. If these signals would be detected, 
these might act as a smoking gun to the bouncing cos- 



mology. 
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